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ABSTRACT 



Aims. To understand the formation of a magnetically dominated molecular cloud out of an atomic cloud. 

Methods. A thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock 
propagating through the hot medium interacts with the cloud. We follow the dynamical evolution of the cloud with a time-dependent 
axisymmetric magnetohydrodynamic code. 

Results. As a fast-mode shock propagates through the cloud, the gas behind it becomes thermally unstable. The ft value of the gas also 
becomes much smaller than the initial value of order unity. These conditions are ideal for magnetohydrodynamic waves to produce 
high-density clumps embedded in a rarefied warm medium. A slow-mode shock follows the fast-mode shock. Behind this shock a 
dense shell forms, which subsequently fragments. This is a primary region for the formation of massive stars. Our simulations show 
that only weak and moderate-strength shocks can form cold clouds which have properties typical of giant molecular clouds. 
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1. Introduction 

Observations of giant molecular clouds (GMCs) show that they 
are magnetically dominated (e.g. Crutcher 1 1999b with values 
of the ratio [5 of thermal gas pressure to magnetic pressure 
of roughly 0.04. Furthermore molecular emission lines show a 
non-thermal broadening component comparable to the Alfven 
speed (Arons & Max 1 1975V However, on scales larger than the 
distances between molecular clouds, the thermal and magnetic 
pressure are about equal. This means that the formation mecha- 
nism of the GMCs must reduce the value of (3 considerably. 

While older stellar associations are devoid of molecular gas, 
molecular clouds in the Solar neighbourhood that do not harbour 
any young stars, are rare. This suggests that the time lag between 
the formation of the clouds and the stars must be short. GMCs 
are thus likely to be transient and dynamically evolving struc- 
tures (e.g. Hartmann, Ballesteros-Paredes & Bergin |2001| >. They 
can be formed by compressive motions of gravitational and/or 
turbulent origin (Ballesteros-Paredes et al. 120061 1. It is thought 
that the formation of GMCs on scales of many tens to hundreds 
of parsecs is regulated by the ram pressure from supersonic flows 
such as supernova remnants, superbubbles and winds of massive 
stars. 

Several groups (e.g. Klein et al. 119941 Mac Low et al. 1 19941 
Fragile et al. 120051 1 studied the interaction of a strong shock with 
a cloud. However, these studies do not concentrate on the for- 
mation of molecular clouds. Instead, they follow the destruction 
of the cloud and its effect on the properties of the interstellar 
medium (ISM). Lim, Falle & Hartquist j2005l hereafter LFH05) 
examined the formation of molecular clouds by the response of a 
warm atomic cloud to an increase in the pressure of the surround- 
ing medium. They implicitly assumed that this pressure increase 
is caused by either flow convergences or weak shocks. Their re- 



sults show that a process relying on pressure-driven compres- 
sion and radiative cooling leads to the formation of highly mag- 
netically dominated regions within a molecular cloud consistent 
with observations of GMCs (Audit & Hennebelle 120051 1. 

In this paper we study the transition of a cloud, triggered by a 
weak shock, from the warm atomic phase to the cold molecular 
one. In Sect. [2] we specify the initial conditions of our model 
and summarise the computational details. Then, we describe the 
dynamical evolution of the cloud (Sect. [3) and present our results 
(Sects. |4] and |5J. Finally, we discuss these results and give our 
conclusions in Sect. [6] 



2. The model 

2. 1 . Initial conditions 

Like LFH05, we assume that a cloud is initially in the ther- 
mally stable warm phase and in pressure equilibrium with the 
surrounding hot ionised gas. The hot gas itself is not in ther- 
mal equilibrium as it is continuously reheated by e.g. supernova 
explosions and superbubbles. Therefore, the net cooling time is 
long. For simplicity, we assume that the gas behaves adiabati- 
cally. Within the cloud, however, heating and radiative cooling 
are important. 

To describe the thermal behaviour of the gas in the cloud, we 
adopt a simplified piece-wise power-law cooling function and 
a constant background heating rate consistent with the standard 
equilibrium phase diagram (i.e. thermal pressure versus density) 
of Wolfire et al. d 1995I ). The net heating rate per unit volume is 
given by 
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H = p[0.015 - pA(T)] erg cm -3 s" 1 , 
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Table 1. Shock Mach Number for the numerical models 



3. Dynamical evolution 



Model 


Shock Mach Number 


A 


2.5 


B 


1.5 


C 


5.0 


D" 


2.5 


E° 


2.5 



" Medium and large cloud models. 
More details are given in Sect.|U 



where 



A(T) = 



3.564 x 10 16 r 212 
9.1 x io 18 r 
1.14 x io 20 r - 56 

1.924 x 10 8 r 3 67 
1.362 x 10 29 r-°- 5 



0K<T< 141 K 
141 K < T < 313 K 
313K<T<6102K 
6102 K < T < 10 5 K 
T > 10 5 K 



(Sanchez-Salcedo et al. 120021 LFH05). For this thermal model, 
the warm phase does not exist for a number density and thermal 
pressure higher than 0.5 cirT 3 and 3051 A;, respectively. Here, k 
is the Boltzmann constant. 

In our calculations, a cloud with initial radius R c i = 200 pc 
has a number density of 0.45 cm -3 . As the cloud is in the warm 
phase, its pressure and temperature, respectively, correspond to 
2825 k and 6277 K. The surrounding hot gas has a number den- 
sity of 0.01 cirT 3 . (Note that these values are the same as in 
LFH05.) To ensure pressure equilibrium between the cloud and 
its surroundings, the temperature of the external medium is set to 
282 500 K. The initial magnetic field is uniform and in the axial 
direction. It has a value such that the ratio of thermal pressure to 
magnetic pressure is unity inside the cloud. 

While LFH05 examined a cloud subjected to a sudden in- 
crease in the external pressure, we consider a steady, planar 
shock hitting the quiescent cloud. This is known as the small- 
cloud approximation (see Klein et al. I19941 >. This intercloud 
shock propagates along the magnetic field lines through the 
hot ionised gas at a speed v ext (quantified by its Mach num- 
ber M = v ext /c s , where c s is the nonmagnetic adiabatic sound 
speed). Table[T]lists the shock Mach number for each model. We 
include models D and E in order to examine how deviations from 
the small-cloud approximation affect the results. 

2.2. Computational details 

All of our calculations were done with an adaptive mesh 
refinement code solving the ideal magnetohydrodynamic 
(MHD) equations for an axisymmetric geometry (Falle & 
Giddings 1993)- The code uses an hierarchy of n levels such 
that the mesh spacing is Ax/2", where Ax is the spacing of the 
coarsest grid. The code refines the grid only where higher resolu- 
tion is necessary. The basic algorithm is a second-order Godunov 
scheme (Falle 119911 ) with a linear Riemann solver and a diver- 
gence cleaning algorithm (Dedner et al. 120021) . 

The computational domain is < r/R c i < 2 and -2 < 
z/Rd S 3. The shock propagates in the negative z-direction and 
we impose fixed boundary conditions, i.e. the shock's down- 
stream values, on the z/R c i = 3 plane. The other boundary con- 
ditions are symmetry on the axis and free flow everywhere else. 
There are six levels of grids with the finest being 1280 x 3200 
which gives a mesh spacing of 0.3 1 pc for our models. Although 
this resolution is high enough to show the onset of the forma- 
tion of dense clumps in the interstellar cloud, it is insufficient to 
follow the complete evolution of the cloud. 



Although the dynamical evolution of the cloud depends on the 
initial shock Mach number, it is similar for all our models. The 
results differ only in detail. Therefore, it is useful to describe 
the evolution before presenting the numerical results in Sects. [4] 
and|5] 

3. 1 . The flow outside the cloud 

When the intercloud shock first hits the cloud, it transmits a 
shock into the cloud and reflects a shock in the previously 
shocked intercloud medium. The reflected shock eventually 
forms a bow shock (or bow wave) in front of the cloud, reduc- 
ing the flow speed of the incoming gas. After a shock-crossing 
time t sc = 2R c i/v ext , the intercloud shock has swept around the 
cloud and converged on the axis behind it. Since the curved, con- 
verging part of the intercloud shock is moving at a lower speed 
than the straight part, swept-up material is raised to a lower pres- 
sure. The transmitted shock is, therefore, weaker at the sides 
of the cloud than on the front and back. As a consequence, the 
shock mainly reduces the cloud size in the direction parallel to 
the distant upstream flow velocity making it oblate. Eventually, 
the cloud ends up as a thin disk. The ram pressure of the external 
flow at the front of the cloud even accelerates this process. 

The transmitted shock has a characteristic speed u,„, = 
Vext/x 1 where y is the density ratio of cloud/intercloud gas. 
This is considerably lower than the propagation speed of the 
intercloud shock. Hence, a velocity shear layer forms at the 
edge of the cloud. This slip surface rolls up to produce a vortex 
ring sweeping material away from the cloud and is also subject 
to Kelvin-Helmholtz instability (Klein et al. 119941 ). However, 
linear instability analysis (e.g. Chandrasekhar 119611 ) predicts 
that a magnetic field strongly suppresses the onset of Kelvin- 
Helmholtz instability. This is also the case for Rayleigh-Taylor 
instability that is excited at the cloud surface due to the acceler- 
ation of the cloud by the rarefied surrounding gas. Although the 
inhibition of these instabilities limits the disintegration of the 
cloud (as seen in hydrodynamic simulations), some fragments 
still get torn off which survive as small-scale coherent structures 
(MacLowetal. [T994j i. 

Because the external flow is diverted around the cloud, the 
magnetic field structure changes significantly. When the shock 
engulfs the cloud, the intercloud gas near the most upstream 
point of contact with the cloud is swept around the cloud and 
carries the magnetic field with it. Since the magnetic field lines 
are anchored in the cloud, they are stretched into a self-reversed 
layer (Mac Low et al. 119941 ). Eventually, these lines double over 
and form a region of strong reversed field. This region may be 
subject to magnetic reconnection. 

3.2. The flow inside the cloud 

A fast-mode shock and a trailing slow-mode shock are driven 
into the cloud. The fast-mode shock compresses the gas so that 
it is in the thermally unstable temperature range. As the time 
scale for the radiative cooling is much shorter than the time for 
the shock to be reflected at the axis, given by the cloud-crushing 
time f cc = Rci/Vim, the gas loses a significant fraction of its in- 
ternal energy during the compression. This initiates the forma- 
tion of cold dense regions embedded in a diffuse warm medium. 
These clumps continue to cool and accumulate material until the 
cooling balances the heating. The clumps end up in approximate 
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Fig. 1. Grey-scale plots of the number density (logarithmic, a) and ft (linear, b) at the times shown for model A. For the number density the range is -2 
The size of the computational box shown is < r/R c / < 1 and -1 < z/R c i < 1. The post-shock flow is in the negative z-direction. 



2 and for /Jit is 0- 1. 



4 



Van Loo et al.: Shock-triggered formation of magnetically-dominated clouds 




Van Loo et al.: Shock- triggered formation of magnetically-dominated clouds 



5 



pressure equilibrium with their surroundings. The number 
density of a clump is then in the range 500 - 1 000 crrT 3 which 
is typical for translucent clumps in GMCs (e.g. Williams et 
al. [19951 

Furthermore, as the magnetic pressure increases across the 
fast-mode shock, the value of /3 falls below unity. This provides 
the ideal conditions for the generation, by MHD waves, of dense 
cores within these clumps (e.g. Falle & Hartquist 2002;; Van Loo, 
Falle & Hartquist I20061 l. especially as thermal instability pro- 
duces a supersonic velocity dispersion within the cloud (Koyama 
& Inutsuka I20021 >. However, we do not have sufficient spatial 
resolution to follow the complete clump and dense core forma- 
tion process. Moreover, self-gravity, which is not included in our 
models, becomes dynamically important once some dense cores 
are compressed to their measured number densities. 

The slow-mode shock moves much more slowly than the 
fast-mode shock. The magnetic pressure decreases behind this 
shock, and consequently there is nothing to stop the gas from 
compressing while it cools. Hence, a shell of cold, dense gas 
forms close to the edge of the cloud. The dense shell first de- 
velops on the upstream parts of the cloud, even before the shock 
has interacted with the rest of the cloud. Furthermore, the veloc- 
ity shear at the surface of the cloud leads to a broader dense shell 
near the stagnation point (see e.g. Fig. [TJ as the gas compressed 
by the slow-mode shock is swept around the cloud. The thick- 
ness of this shell thus depends on the level of shear at the edge, 
i.e. the shock Mach number. 

The interaction of the shock with the cloud produces dense 
regions in the cloud, i.e. a shell at the edge of the cloud and 
several individual clumps within the cloud. The shell contains a 
large fraction of the cloud's total mass. It is also unstable due to 
thermal instability and dynamical instabilities (such as Kelvin- 
Helmholtz and Rayleigh-Taylor), so that it fragments into cold, 
massive filaments. Our simulations show this fragmentation at 
later times. This means that the shell may well be the precursor 
of massive star formation. The clumps within the cloud, may 
form only low-mass stars as they are magnetically dominated. 
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Fig. 2. The fraction of the cloud volume with /3 < 0.1 as a func- 
tion of time (in Myr) for Model A (solid), Model B (dashed) and 
Model C (dashed-dotted). 



,,0 T 

0.9- 
0.8- 
0.7- 

[Z 

O 0.6- 




2 4 6 8 10 12 14 16 18 20 

t (Myr) 



4. Numerical results 

Figure Q] shows grey-scale plots of the number density and /3 for 
model A. We can clearly see that, as the fast-mode shock moves 
through the cloud, the swept-up gas becomes magnetically dom- 
inated, i.e. (3 drops below unity. This is because the thermal gas 
pressure decreases due to the transition to the thermally unsta- 
ble phase, while the magnetic pressure increases behind the fast- 
mode shock. It can then be expected that, by the time the shock 
reflects at the symmetry axis, a significant fraction of the vol- 
ume of the cloud has a low-/3 value. Figure [2] shows the fraction 
of cloud volume with fi < 0.1 as a function of time for each of 
three models. At around f„ , the volume fraction of the model 
A cloud for which /3 < 0.1 is about ~ 0.9. Actually, most of 
the cloud is in a low-/? regime for a significant period of time 
(~ 3 Myr). Within this time range, the cloud has a characteris- 
tic mean number density and radius of w 20 crrT 3 and * 50 pc. 
Furthermore, we find that the velocity dispersion of the cloud is 
highly supersonic. These are all typical properties of a GMC. 

Although a small fraction of the volume of cloud is not mag- 
netically dominated, its corresponding mass fraction is consid- 
erable (see Fig.0. About half of the total mass is in the bound- 
ary layer of the cloud. Here, the gas is in the cold phase at a 
high thermal pressure so that /? « 1. While the conditions in 
the interior of the cloud are ideal for the generation of dense 
clumps and embedded cores by MHD waves, this is not true for 



Fig. 3. The fraction of the cloud mass with (3 < 0.1 as a function 
of time (in Myr) for Model A (solid), Model B (dashed) and 
Model C (dashed-dotted). 



the dense shell. It is, however, expected that this shell breaks 
up into dense fragments (Koyama & Inutsuka |2000| >. This frag- 
mentation is observed in our simulations. The densest fragments 
have number densities of ~ 10 3 cirT 3 and sizes of ~ 1 pc. With 
a temperature of about 30 K for these fragments, we derive a 
Jeans mass of ~ 350 Mq which is similar to the inferred mass. 
(We have assumed a near spherical shape for the dense frag- 
ments.) Therefore, the boundary layer may be a site of massive 
star formation. 

For the other models, we find similar results, except that 
stronger shocks (Model C) lead to a faster evolution (see Fig.[2j. 
This is because the propagation speed, u,„ f , of the shock in- 
creases with Mach number, reducing the dynamical time-scale 
t cc . (The cloud crushing time is t cc k 4 Myr for Model C, while 
t cc « 8.5 Myr for Model A.) Magnetically-dominated regions 
thus develop on shorter time-scales. However, the compression 
of the cloud into a disk also proceeds faster due to a higher ram 
pressure exerted by the external flow (i.e. the ram pressure is 
proportional to M 2 ). Our simulations show that, for Model C, the 
time-scale for flattening is ~ lt cc , while it increases to ~ 1 .5-2t cc 
for weaker shocks (Model A and B). Consequently, the lifetime 
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of a molecular cloud produced by a strong shock has an upper 
limit of a few Myr before it is flattened. The formation of H2 
and CO from atomic gas requires a similar timescale (Bergin et 
al. 120041 Glover & Mac Low |2007l >. Thus, as the dynamical evo- 
lution is too rapid for hydrogen to be completely converted to H2, 
strong shocks do not produce giant molecular clouds. However, 
molecules can still form in the flattened cloud. As this thin disk 
is being shredded by instabilities (MacLow et al. |1994| l, this re- 
sults in the formation of small molecular clouds with sizes of a 
few pc. 

The time-scale for molecular cloud formation by a weak 
shock (Model B) is compatible with the expected time-scales for 
the formation of molecules (the cloud flattens on a time-scale of 
~ 25 Myr). However, in this case the fast-mode shock propagat- 
ing through the cloud does not compress the gas sufficiently to 
induce fast radiative cooling. As the cooling time of the gas is 
longer, it is difficult to produce magnetically-dominated regions 
within the cloud. Figure|2]shows that less than 10 % of the cloud 
has p < 0.1 during its entire evolution. Note, however, that we 
have not included any feedback effect of star formation. It is be- 
lieved to take place in the dense boundary layer of the cloud. We 
can then expect proto-stellar outflows to reduce the value of /3 
locally, thus, forming magnetically dominated regions. 



5. Medium and large clouds 

In the previous section we showed the results for a shock, being 
driven by a constant distant upstream pressure, interacting with 
a cloud; i.e. the so-called small-cloud approximation. However, 
the pressure behind a wind-blown bubble or a supernova blast 
wave can change significantly on a time-scale shorter than t cc , 
or even t sc . By comparing these time-scales to the time-scale, 
t p , of pressure variations we can make a further distinction be- 
tween medium-sized and large clouds (Klein et al. [1994)l. For 
medium-sized clouds, the pressure does not vary much as the 
shock sweeps around the cloud, but it does change significantly 
before the cloud is crushed, i.e. t cc > t p > t sc . For large clouds, 
t p < t sc , and the pressure thus changes while the shock sweeps 
around the cloud. 

To study the effect of pressure variation behind the shock, 
we assumed that a shock is accelerated exponentially from M = 
1.1 to a higher value of M (given in Table Q] for models D and 
E) after which it accelerates no more. In our calculations, we 
assumed that the period of constant shock strength is ~ 1 Myr 
for a large cloud (Model E) and ~ 5 Myr for a medium-sized 
cloud (Model D). Our model is a simplified description of a 
blast wave. However, we prefer this to the implementation of the 
Sedov-Taylor blast-wave solution, as we are only interested in 
the qualitative differences from the small-cloud approximation. 

Figure [4] shows that, for medium-sized clouds, the results 
do not change significantly from the small-cloud approximation. 
Large clouds, however, evolve at a lower rate towards a low-/? 
regime. This can be readily explained by the fact that the pres- 
sure changes as the shock sweeps around the cloud so that the 
shock strength on the upstream and downstream sides is quite 
different. The fast-mode shock therefore does not produce the 
same compression of the gas as in the small-cloud approxima- 
tion. This results in a longer cooling time and a delay in the for- 
mation of a low-/? cloud. Yet, most of the cloud in the large-cloud 
approximation still becomes magnetically dominated. 

The fact that the pressure behind the shock decreases on 
time-scales shorter than t sc has the additional advantage that 
the cloud is exposed to a lower ram pressure. Consequently, the 
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Fig. 4. Similar to Fig. |2 but now for small (solid), medium 
(dashed) and large (dotted) cloud approximation. 




Fig. 5. Logarithmic grey-scale plot of the number density for 
Model A (left) and Model E (right) at 10 Myr. The range for 
the number density is -2-2. 



cloud does not flatten into a disk as fast as in the small-cloud ap- 
proximation (see Pig . ■ This is relevant as it shows that shocks, 
stronger than the small-cloud approximation would indicate, can 
produce clouds resembling GMCs. 

6. Conclusions and discussion 

In the present paper we studied the interaction of a shock with 
an initially warm, thermally stable cloud of uniform density. 
As the shock sweeps around the cloud, a transmitted fast-mode 
shock propagates through the cloud producing magnetically- 
dominated regions behind it. This provides the ideal condi- 
tions for MHD waves to produce high-density clumps within the 
cloud. These clumps may eventually form low-mass stars. 

A slow-mode shock follows the fast-mode shock leading 
to a high-density layer near the boundary of the cloud. As the 
shock interacts first with the front of the cloud, the dense layer 
forms there first. The dense layer is subject to instabilities and 
fragments into massive clumps. Our analysis shows that these 
clumps are close to gravitationally bound and, thus, are possible 
precursors of massive stars. This agrees well with the findings 
that OB stars form almost always near the boundary of molecular 
clouds and even prefer to form at one particular edge of the cloud 
(e.g. Blaauw[T9g2l El megree n & Lada [1977l Israel [T978l Gatley 
et al. 119791 Fich et al. [1982). Furthermore, the cometary cloud 
structure seen in our simulations, i.e. a massive head and long- 
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Galactic longitude 



Fig. 6. The W3 GMC in the J=l -0 tran sition of 12 CO with a 50 
arcsecond resolution (Bretherton[2003). High-mass star-forming 
regions are found in the HDL situated in the North-East of the 
cloud. 



spread tail, seems characteristic for clouds harbouring cluster- 
forming cores (Tachihara et al. l2002l >. 

In our simulations we imposed a sharp boundary between 
the cloud and the surrounding medium. However, real interstel- 
lar clouds have an internal density gradient, e.g. the southwest 
cloud in the Cygnus Loop supernova remnant shows a grad- 
ual drop in density near the boundary (Patnaude et al. 120021 ). 
Hydrodynamical simulations of the effect of smooth boundaries 
show that the initial density pattern remains imprinted on the 
shocked cloud (Nakamura et al. 2006). Preliminary results of 
MHD simulations show a similar effect, but we still find that a 
thin dense layer arises near the boundary of the cloud. The effect 
of a density gradient and the presence of condensations within 
the cloud will be investigated in detail in a subsequent paper. 

A likely example of a shock-cloud interaction is found in the 
W3 GMC. Figure shows the W3 GMC as traced b y 12 CO J=l- 
emission at 50 arcseconds resolution (Bretherton 120031 . The 
cloud is ~ 60 pc in size with a total gas mass of nearly 4 x 10 5 
M , of which ~ 40% is in a boundary layer referred to as the 
high-density layer (HDL) (Lada et al. 119781 ). The mean number 
density in the cloud is around 40 cirT 3 implying a very clumpy 
medium in the cloud. A number of high-mass star-forming re- 
gions are found in the HDL which runs parallel to the edge of 
the W4 H ii region and is likely to be formed by the expansion 
of the H ii region and the stellar winds of the W4 OB associa- 
tion. Recently, Moore et al. (120071 ) surveyed W3 and found that 
the fraction of the total cloud mass in dense clumps is as high as 
0.26 in the HDL and only 0.05 in the diffuse part of the cloud. 
These values are consistent with those obtained in our simula- 
tions. 

Our simulations support the idea that molecular clouds are 
transient objects. In this view, star formation happens quickly 
and the parental GMCs are short-lived with observations sug- 
gesting lifetimes of the order of 10 - 20 Myr (see e.g. Blitz 
et al. 120071 ). A lower limit of a few Myr is set by the forma- 
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tion of H2 and CO from the atomic gas (Bergin et al. 120041 ). 
Although the lifetimes of the cloud in our simulations can only 
give a rough guide because the cloud compression into a disk 
is fastest when the external flow is along the magnetic field 
lines, they agree well with observed values. For strong shocks 
(M > 5), however, the evolution seems to be too fast to produce 
giant molecular clouds with morphologies as those observed. 
Nevertheless, strong shocks can still produce molecular clouds 
as the formation of molecules continues within the small coher- 
ent structures that fragment off the initial cloud. 

Our results would also seem to rule out very weak shocks as 
agents for molecular-cloud formation since they do not produce 
clumped clouds. A deficiency of our model is that it does not 
include any feedback from star formation such as proto-stellar 
winds and jets, and expanding H 11 regions. However, follow- 
ing the sequential star formation picture of Elmegreen & Lada 
(11 9771 ), we argue that these processes occur in the boundary layer 
of the cloud and inject sufficient energy into the cloud to produce 
a clumped medium. We, therefore, conclude that only weak to 
moderate-strength shocks are the likely triggers for the forma- 
tion of GMCs in the ISM. 
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